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ABSTRACT 

Galaxy clusters can be detected as surface brightness enhancements in smoothed optical surveys. This method 
does not require individual galaxies to be identifiable, and enables clusters to be detected out to surprisingly high 
redshifts, as recently demonstrated by the Las Campanas Distant Cluster Survey (LCDCS). Here, we investigate 
redshift limits for cluster detection in the Sloan Digital Sky Survey (SDSS). Calibrating assumptions about the 
surface brightness profile, the mass-to-light ratio, and the spectral energy distribution of galaxy clusters using 
available observational data, we show that it should be possible to detect galaxy groups out to redshifts of ^ 0.5, 
and massive galaxy clusters out to redshifts of ^ 1.2 in summed / + i' + z! SDSS data. Redshift estimates can 
be derived from the SDSS magnitudes of brightest cluster members out to redshifts near unity. Over the area 
of sky it covers, SDSS should find > 98% of the clusters detectable by the Planck satellite through the thermal 
Sunyaev-Zel'dovich effect. The few Planck clusters not detected in SDSS will almost all be at z > 1.2. 



1. Introduction 



I 

o 



> 

Galaxy clusters are not only interesting for studying galaxy evo- 
lution; they have also become an important tool for understand- 
ing the growth of cosmic structure and the cosmological frame- 
T— I work in which it occurs. The construction of observed galaxy 
T— I cluster samples out to redshifts unity and beyond is thus an im- 
portant goal of current observational cosmology. Traditionally, 
galaxy clusters have been identified at optical wavelengths as 
regions on the sky where the number density of galaxies suffi- 
- ciently exceeds its mean. This definition requires of course that 
the cluster galaxies be individually detectable, setting an upper 
imit to redshifts at which galaxy clusters can be found. 

Dalcanton (1996) proposed that clusters could be detected as 
regions on the sky where the surface brightness exceeds the av- 
erage sky brightness. Her suggested procedure consists of re- 
moving galaxies from carefully flat-fielded images, smoothing 
the residual image with a kernel whose width should approxi- 
mately match the angular extent of galaxy clusters at intermedi- 
ate redshifts, and searching for peaks in the smoothed surface- 
brightness distribution which sufficiently exceed the noise level 
of the smoothed sky background. 

Gonzalez et al. (2001) successfully applied this technique to 
the Las Campanas Distant Cluster Survey (LCDCS) data taken 
with the Las Campanas Great Circle Camera (Zaritsky, Schect- 
man & Bredthauer 1996) and constructed a catalog of 1073 
groups and clusters. Estimated redshift limits of the catalog 
range from ^ 0.3 for groups to ^ 0.8 for massive galaxy clus- 
ters. 

Being intrinsically highly uniform by construction, drift-scan 
surveys like the LCDCS provide an optimal data basis for the 
application of this cluster-detection technique. In this paper, 
we estimate the redshift limits expected for the largest ongoing 
drift-scan survey, the Sloan Digital Sky Survey (SDSS; York et 
al. 2000). Compared to the LCDCS, the SDSS uses a mirror area 
larger by a factor of 6.25, an exposure time shorter by a factor of 
3.6, and a system of broad-band filters rather than a single, very 
broad filter. On the whole, the flux limit of the SDSS is expected 
to be roughly a factor of two below that of the LCDCS, which 
should allow the detection of massive clusters out to redshifts 
beyond unity. We study this expectation in detail in this paper 
We note that cluster catalogues constructed from the SDSS by 
more traditional techniques (i.e. based on galaxy catalogues) are 



already expected to be complete to z ^ 0.4 (Kim et al. 2001). 

Sect. ^ describes our assumptions. Results are presented in 
Sect. ^ and we summarise our conclusions in Sect. ^ 



2. Assumptions 

Throughout, we adopt a flat, low-density cosmological model 
with matter density parameter Q.q = 0.3, cosmological constant 
corresponding to density parameter £2a = 0.7, and a Hubble con- 
stant //o = 70kms^^ Mpc^'. The CDM power spectrum is taken 
to have primordial spectral slope « = 1 and is normalised using 
08 — 0.93 so as to reproduce the local abundance of galaxy clus- 
ters. 

2.1. The SDSS photometric system 

The SDSS photometric system uses five wide, almost non- 
overlapping bands covering the wavelength range between 
3000 A and 11 000 A. The recently re-measured filter response 
functions S^ are shown in Fig. |l] (M. Strauss, private commu- 
nication). The effective wavelengths A,eff, effective band widths 
AXeff, and flux sensitivity quantities 



Q = Jd{\nv)Sy 



(1) 



of the five bands are summarised in Tab. |l]. 

The system is described in detail by Fukugita et al. (1996). 
Its zero points are placed on the spectrophotometric AB mag- 
nitude system, so that magnitudes can directly be converted to 
fluxes in physical units. Given the spectral flux in units of 
ergcm^^s^' Hz^^, the broad-band AZ? magnitude is defined as 



fd(lnv)FySv 
-2.51ogio -^ ^ Q ' -48.60, 



(2) 



with Q defined in Eq. ([I]). 

In order to assess the noise level of the sky background, we 
need to relate the incoming flux to the detector counts. The num- 
ber of photoelectrons released by the incoming spectral flux Fv 
is given by 

-1 



Ne=Ath-' / d(lnv)/v5- 



(3) 
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Fig. 1. — Recently re-measured SDSS filter response functions. From 
left to right: u', g', r', i', and z'- 
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21.3 
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Table 1. — Effective wavelengths A,eff, effective band widths 
AA-eff, and flux sensitivity quantities Q of the five SDSS bands. 
Also listed is the sky brightness in the five SDSS photometric 
bands at the SDSS telescope site at Apache Point, and corre- 
sponding counts. The last row gives the 5-0 Poisson fluctuation 
level of the sky background, /j^ky and A/j^^y are given in magni- 
tudes per square arc second. 



where h is Planck's constant. The survey telescope area is A = 
4.9 X lO'^cm^, the effective exposure time per band is f = 54.1 s. 

The main sources of noise are the count fluctuations of the 
sky brightness and the brightness fluctuations in galaxies just be- 
low the SDSS detection limit. For the sky brightness /jj^y at the 
SDSS telescope site at Apache Point, we use values measured 
by M. Richmond (see Tab. nl given in magnitudes per square 
arc second). Given the area A^of the SDSS survey telescope, the 
effective exposure time t, and the SDSS response functions, the 
brightness of the night sky can be converted to counts, which are 
listed for each band in Tab. |l} The Poisson fluctuation level of 
the night sky per square arc second is on the order of a few per 
cent. The 5-0 Poisson fluctuation in magnitudes per square arc 
second is given in the last row of Tab. ^ 

The noise contributed by undetected galaxies is small com- 
pared to the shot noise from the sky brightness. Because of the 
flatness of their number-count function, most of the light from 
undetected galaxies comes from galaxies just below the detec- 
tion limit. Suppose the typical solid angle covered by such a 
galaxy is A, and their number density per solid angle is tig^i. Let 
the photon count produced per solid angle of the night sky be 
«e,sky, and assume that the detection threshold is v times the sky 
brightness fluctuations. The number of sky photons per galaxy 

is then (we.sky^), and the detection threshold is v(ne,sky^)^''^- 
The surface brightness from galaxies at the detection threshold 



is then 

Sgal = «galV(«e,sky^)'''^ , 

and its Poisson fluctuation in the solid angle 5£2 is 

A5gal = (ngal5n)'/2v(„e,skyA)'/2. 

The Poisson fluctuation in the sky brightness is 

A5sky = («e.sky8ii)'''^ • 



(4) 



(5) 



(6) 



This implies that the Poisson fluctuations in the flux from unde- 
tected galaxies are small compared to the sky brightness fluctu- 
ations, since the ratio 



A5gal 
A5sky 



= v(ngaiA) 



1/2 



(7) 



scales with the square root of the sky fraction covered by galax- 
ies, which is on the per cent level. 

The variance in the flux from undetected galaxies exceeds the 
Poisson expectation because the galaxies are correlated. The 
contribution from the correlation is, however, smaller by ap- 
proximately a factor of ten than the Poisson fluctuations at the 
apparent magnitude levels of relevance (e.g. Brainerd, Small & 
Mould 1995), so that the total fluctuations expected from un- 
detected galaxies are negligible compared to the sky brightness 
fluctuations. 

2.2. Cluster spectra 

We describe the galaxy cluster spectrum as a weighted su- 
perposition of an early- and a late-type galaxy spectrum, and 
F^, respectively. The combined cluster spectrum at redshift z is 
then written as 



(8) 



The early-type fraction of the cluster population, fsiz), de- 
creases with redshift. The results of van Dokkum et al. (2000) 
suggest choosing a linear decrease with redshift. 



/£(z) =0.8-0.4z. 



(9) 



The mass-to-light ratio of early-type galaxies in clusters de- 
creases with increasing redshift. For the mass-to-light ratio 
in the B band, van Dokkum et al. (1998) find Alog(M/LB) - 
— 0.4z. Adopting this relation, we assume that early-type galax- 
ies brighten with redshift as 



Ie{z) = 10' 



OAz 



(10) 



The late-type fraction increases with redshift, but we assume 
the mass-to-light ratio of late-type galaxies to be constant. In 
order to account for the fact that late-type galaxies are typi- 
cally fainter than early-type galaxies, we introduce a redshift- 
independent factor fi < I, which, based on the data of van 
Dokkum et al. (2000), we set to 



fi 



(11) 
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The spectra F^ ' themselves are redshifted, but otherwise as- 
sumed to be constant. They were kindly provided by S. Chariot 
(private communication). 
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Table 2. — Mass-to-light ratios in solar units of a galaxy cluster 
at redshift zero in the five SDSS bands. The mass-to-light ratio 
is normalised to 250 in the Johnson B band. For comparison, the 
mass-to-light ratios in the conventional Johnson bands are also 
shown. Due to the predominantly red galaxy population, the 
mass-to-light ratio decreases towards longer-wavelength bands. 



The mass-to-light ratio of galaxy clusters is typically mea- 
sured to be on the order of M /Lb = ISQMq/Lq in the B band 
(e.g. Bahcall, Lubin & Dorman 1995; Carlberg et al. 1996; Carl- 
berg, Yee & Ellingson 1997a; Gonzalez et al. 2000), however 
with a scatter of ^ 20%. We normalise the combined cluster 
spectrum such that the B-band mass-to-light ratio at redshift zero 
is 250 in solar units. Specifically, let /\p be the solar spectrum, 
then the amplitude of the cluster spectrum is chosen such that 

J d(lnv)if 5v.B = J d(lnv)F®5v.fl (12) 

is satisfied. The mass-to-ratios in all five SDSS bands and in the 
conventional ]ohnson-U BVRI bands are given in Tab. |[ 

The predominantly early-type galaxy population causes 
the mass-to-light ratio to decrease towards longer-wavelength 
bands. This trend agrees well with observations; for instance, 
the /-band mass-to-light ratio of Abell 1651 was measured to be 
- 160 (Gonzalez et al. 2000). 

Figure ^ shows arbitrarily normaUsed cluster spectra at red- 
shifts zero and unity. 
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Fig. 2. — Composite cluster spectra as functions of wavelength X in 

A, at redshifts zero (solid curve) and unity (dotted curve). Both spectra 
are arbitrarily normalised. 

2.3. Cluster k-corrections and colours 

Integrating the synthetic cluster spectrum with the SDSS filter 
curves as in Eq. (0), we can now compute ^-corrections and clus- 
ter colours for the SDSS photometric system. Figure || displays 
the cluster ^-corrections in the five SDSS bands as functions of 
cluster redshift. 
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Fig. 3. — Cluster A;-corrections as a functions of redshift for M = 5 x 
IO'^'Mq//! in the five SDSS bands (solid curve: u , dotted curve: e , 
short-dashed curve: /, long-dashed curve: dash-dotted curve: z') 



The ^-correction in the bluest SDSS band (m') grows to ~ 
2.4 mag. up to z ^ 0.5 and then levels off. This is because at that 
redshift essentially all of the cluster flux has been shifted red- 
ward of the m' filter curve, and the remaining short-wavelength 
part of the cluster spectrum is almost flat. The k corrections in 
the other SDSS bands grow monotonically with redshift. Their 
amplitudes decrease as the filters get redder 

Figure ^ shows the four different cluster colours iu'-g'), ig'- 
r'), {r'-i') and (i'-z!) as a function of cluster redshift. 
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Fig. 4. — Cluster colours as functions of redshift, for cluster mass M = 
5 X 10'^M,;,//j. Solid curve: (u'-g'); dotted curve: (g'-r'), short-dashed 
curve: {r'-i'); long-dashed curve: (i'-z')- 

On the whole, clusters become redder as their redshift in- 
creases, with the exception of the bluest colour shown (u'-g'), 
which starts red at low redshift and then decreases beyond red- 
shift z ^ 0.4 as the 3000 A break shifts through the g' band. 

2.4. Cluster surface brightness 

We assume that the average light distribution in galaxy clusters 
follows mass, and that the mass distribution is described by the 
familiar NFW density profile (Navarro, Frenk & White 1996, 
1997). Data from the CNOC survey shows that this gives a very 



good model for the mean observed luminosity profile of rich 
clusters (Carlberg et al. 1997b). The three-dimensional matter 
density is 

where is a characteristic scale radius, Ps is the density scale, 
and x = r/rs. Although two free parameters appear in ([13), there 
is effectively only one parameter to set because the scale radius 
and the density scale are related. We choose the virial mass M as 
the remaining single free parameter and compute rs and Ps from 
it, following the prescription given in Navarro et al. (1997). The 
virial mass is 



M = 47tpsrs 



ln(l 



1- 



(14) 



where the concentration c = rioo/rs is the ratio between the 
virial radius and the scale radius. 

Given the cluster mass, its luminosity L is determined by the 
mass-to-light ratio 

M 

When projected along the line of sight, its surface mass density 
of an NFW halo is E(x) — Ps r^ f(x), with 



fix) 



1 



: arctan 



1 



.T-l 
X+l 



(16) 



(see e.g. Bartelmann 1996). Therefore, the surface brightness 
profile of the cluster is 



Scix) 



1 -1 



ln(l 



1 



(17) 



Sc{x) is the energy radiated by the cluster per unit unit surface 
area per unit time. 

Photon conservation, and the Etherington relation between 
angular-diameter and luminosity distance then imply that the 
flux received from a cluster at redshift z per unit time, frequency, 
detector area and solid angle is 



47i:(l 
D(z)0 



z)^4jtr? 



ln(l 



/ 



(18) 



where 9 is the angular separation from the projected cluster cen- 
tre, D{z) is the angular diameter distance to the cluster, and 
is the cluster spectrum of Eq. (^, i.e. the cluster luminosity per 
unit frequency. 

By means of Eq. (0), Sy{Q) can now be converted into a sur- 
face brightness profile for the cluster in the conventional units 
of magnitudes per square arc second. The dotted curve in Fig. ^ 
shows an example for a cluster of mass 2 x lO'^^M©//! at redshift 
z = 0.75. 



2.5. Smoothing 

The signal-to-noise ratio of a cluster detection can be increased 
by smoothing. We therefore convolve the projected cluster pro- 
file ( p^ ) with a Gaussian filter of width A9. Exploiting the ax- 
ial symmetry both of the cluster profile and of the Gaussian, 



and making use of the convolution theorem in Fourier space, the 
convolution can be transformed into 



5v"(e) = 



Ae2 



exp 



2Ae2 



(19) 



where Io(x) is the modified Bessel function of order zero. It 
can be approximated by lo(x) « (27tx)^'/^ exp(x) for l,in 
which case the convolution simplifies to 



exp 



(9-^) 



21 



2A92 



(20) 



We show in Fig. || the central surface brightness in the r' band 
of a cluster with mass M = 5 X IQ^'^MQ/h at redshift 0.75, 
smoothed with a Gaussian kernel of varying width A9, and the 
5-0 sky brightness fluctuation level. The small inserted panel 
shows the difference between the two curves. 
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Fig. 5. — Central surface brightness in the r' band (in magnitudes per 
square arc second) as a function of smoothing radius for a cluster of 
mass M = 5 X 10^* Mpj/h at redshift z = 0.75 (solid curve) and the 5-0 
sky brightness fluctuations (dotted curve). The small panel shows the 
difference between the curves. 

Figure |] shows that the sky noise falls below the central clus- 
ter surface brightness when the smoothing kernel grows above 
« 0.2', and the height of the cluster centre above the noise 
reaches a broad maximum at A9 w 0.5'. Correspondingly, we 
choose a smoothing kernel width of A9 = 0.5' in the following. 

Figure ^ illustrates the effects of smoothing and sky noise 
fluctuations on the surface-brightness profile of a cluster with 
mass M = 2 X lO'^^M© /h at redshift z — 0.75. The intrinsic pro- 
file is broadened by the Gaussian smoothing kernel, which has 
A9 = 0.5'. While the smoothed profile falls below the 5-0 noise 
fluctuation level if data in the r' band only are used, the cluster 
is well detected if data in the g', r' and ;' bands are summed. 

3. Results 

We can now proceed to compute the upper redshift limit for a 
significant cluster detection. We first consider individual SDSS 
bands. Figure ^ displays the redshift Zmax as a function of clus- 
ter mass at which the central cluster surface brightness in each 
of the five SDSS bands drops below the 5-0 limit of the sky 
background fluctuation, after smoothing with a Gaussian kernel 
of width A9 = 0.5'. 
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Fig. 6. — Surface brightness profiles (in magnitudes per square arc sec- 
ond) for a cluster of mass M = 2 x 10''*M©//i at redshift z = 0.75 as 
a function of cluster-centric distance 6. The assumed intrinsic NFW 
surface brightness profile in the r' band is shown as the dotted curve, 
which is smoothed by a Gaussian kernel of width A6 — 0.5' (dashed 
curve). The solid curve shows the surface brightness summed in the g', 
r' and /' bands. The horizontal curves indicate the 5-c sky brightness 
fluctuation levels in the r' band (dashed) and the summed g', r' and 
bands (solid), respectively. While the smoothed cluster falls below the 
noise level in the r' band, it is well detected in the summed bands. 



Clearly, the detection limit in the u' band is the poorest, reach- 
ing only out to redshift z ^ 0.15 at M = lO'^^M©//!. This is 
a consequence of the combined effect of the relatively high 
background fluctuation level, the low detector efficiency and 
the predominantly red cluster light. The detection limit in- 
creases rapidly with increasing filter wavelength to reach just 
above redshift unity in the /' band for massive clusters with 
M = 10^^ M(v,/h. The relatively low k correction in the z' band 
leads to the comparatively steep increase in the upper redshift 
limit despite the poor efficiency in the redmost SDSS band. 

Figure 1^ indicates that the detection limit can be considerably 
increased by summing data in several bands. Figure ^ quantifies 
this expectation. 

The solid curve shows the detection limit for data summed 
in the two bluest bands (m' and g'), the dotted curve for data 
summed in the two redder bands (/ and /'). The increase be- 
tween the two is substantial; the upper redshift limit approxi- 
mately doubles. One to two tenths in redshift are gained if the 
z' band data are added to the r' and /' bands. Further addition of 
the g' band improves the limit a little more, while the relatively 
low signal-to-noise ratio in the u' band even lowers the detection 
limit if data in all bands are summed. 

These results show that massive clusters with M ^ lO^^Mpj/Zi 
should be detectable as significant surface-brightness enhance- 
ments in the SDSS data out to redshifts of z ^ 1.25 if data 
in the three redmost bands are summed. The redshift limit 
of the cluster detection increases approximately linearly with 
logio(M/i/M0). 

This redshift limit agrees with expectations raised by the Las 
Campanas Distant Cluster Survey (LCDCS; Zaritsky et al. 1997; 
Gonzalez et al. 2001). The LCDCS used a 1-m telescope with 
an effective exposure time of 194 s, while the SDSS has a 2.5-m 
telescope and an effective exposure time of 54.1 s. The product 
of telescope area and exposure time is therefore larger for the 
SDSS by 75%. The LCDCS uses the broad W filter. It covers 
the wavelength range between 4600 A and 7300 A, which over- 



10" ID'* 10' = 

M [Mg/h] 

Fig. 7. — 5-0 detection limits for clusters as a function of mass in the 
five individual SDSS bands. Solid curve: u', dotted curve: g', short- 
dashed curve: r', long-dashed curve: dash-dotted curve: z'. 
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Fig. 8. — 5-CT detection limits for clusters as a function of mass in sev- 
eral combinations of wave bands. Solid curve: (t/ + g'); dotted curve: 
(r' + i'); short-dashed curve: (r' + 1' +z'); long-dashed curve: all bands 
except ii\ dash-dotted curve: all bands. 



laps with the SDSS r' filter and major fractions of the ^ and 
/' filters. The combined width of the + r' + /' +z' filters is 
approximately 4000 A, or ^ 50% larger than the width of the 
W filter. Furthermore, adding z' data helps because of the pre- 
dominantly red colour of galaxy clusters. Adopting the same 
quantum efficiency for the SDSS and LCDCS CCDs, it seems 
reasonable to assume that the combined ^' + r' + /' +z' SDSS 
data lower the LCDCS flux limit by approximately a factor of 
(1.5 X 1.75)'/2w 1.6. 

The LCDCS contains groups out to redshifts of ~ 0.3, and 
massive clusters out to redshifts of ~ 0.8, while we estimate 
respective redshift limits for SDSS to be ^ 0.5 and ~ 1.3. As 
Eq. ( [l8| ) shows, the observed surface brightness 5^ scales with 
redshift roughly as (1 +z)^^. The estimated improvement in 
the flux limit by ^ 1.6 thus leads to an expected increase in the 
redshift limit of 

zLT/^'^ 1.6^/^(1 (21) 
which agrees reasonably well with our direct estimates. We thus 



confirm the speculation by Gonzalez et al. (2001) that the SDSS 
redshift limit for massive clusters could reach ~ 1 .25. 

We note that it should be possible to estimate redshifts for 
a large fraction of the SDSS cluster sample through the mag- 
nitudes of the brightest cluster galaxies (BCGs). They form a 
remarkably homogeneous class of objects (Hoessel, Gunn & 
Thuan 1980; Schneider, Gunn & Hoessel 1983; Postman & 
Lauer 1995) with a very narrow luminosity function. The K- 
band absolute magnitudes of BCGs in the redshift range < 
z < 1 are compatible with no luminosity evolution, and their 
colour evolution is consistent with a passively evolving, old stel- 
lar population (Aragon-Salamanca, Baugh & Kauffmann 1998). 
Assuming a reduced absolute magnitude of 21.1 in the g' band 
(cf. Schneider et al. 1983) and a non-evolving early-type spectral 
energy distribution, we estimate the apparent BCG magnitudes 
shown as functions of redshift in Fig. p[ We checked and con- 
firmed that they are compatible with the /T-band measurements 
by Aragon-Salamanca et al. (1998) 
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Fig. 9. — Estimated magnitudes of brightest cluster galaxies in the g', r' 
and z' SDSS filter bands (dotted, short-dashed, and long-dashed curves 
as indicated), and in the summed filters (solid curve). The dash-dotted 
curve indicates a conservative detection limit of 22.5 magnitudes (ap- 
proximately 5-<J). Brightest cluster galaxies should be detectable in 
summed SDSS data out to redshifts around unity. When identified with 
galaxy cluster candidates, they will increase the reliability of the detec- 
tion and allow photometric redshift estimates. 

If data in the g', r' and /' data are summed, it should be pos- 
sible to identify BCGs out to redshifts near unity, assuming a 
conservative detection limit of 22.5 in the summed bands. Note 
that the estimated r' magnitudes out to redshifts of ^ 0.5 agree 
very well with the r' magnitudes of galaxies with measured red- 
shifts in the SDSS Luminous Red Galaxy Sample (Eisenstein et 
al. 2001). 

Upcoming wide-area surveys in the sub-millimetre regime, 
like that planned with the Planck satellite (Bersanelli et 
al. 1996), will allow the detection of galaxy clusters through 
their peculiar spectral signature caused by the thermal Sunyaev- 
Zel'dovich effect. Taking the projected temperature sensitivity 
of Planck and its characteristic beam width of 5' as a baseline 
and modelling the cluster population with standard assumptions, 
Bartelmann (2001) estimated the upper redshift limit of the clus- 
ter sample expected to be detectable for Planck. Cumulative 
redshift distributions for clusters with M > 5 x lO^^/i"' de- 
tected in SDSS and Planck data are shown in Fig. |o[ 

As explained in Bartelmann (2001), Planck is expected to 
see clusters at any redshift provided the cluster mass exceeds 




Fig. 10. — Cumulative redshift distributions A'(> z) of clusters with 
M > 5 X 1O'^/?"'M0 in summed r'+Z'+z' SDSS data, detected at 2- 
CT and 5-0 significance as indicated in the plot. Also shown is the 
cumulative redshift distribution of clusters detected through the ther- 
mal Sunyaev-Zel'dovich effect at 2-o significance by the Planck satel- 
lite. Shown is the expected number of clusters in the SDSS area. 
The abundance of SDSS clusters lies almost two orders of magnitude 
above that for Planck clusters up to the redshift where clusters with 
M = 5 X 10^^ hr^ Mq:; drop below the detection limit. The SDSS abun- 
dances then cut off exponentially. All Planck clusters in the survey area 
are detected by SDSS out to z --^ 1.1 at 5-CT, and out to z --^ 1.3 at 2-o. 
Cluster catalogues constructed from the SDSS data by more traditional 
techniques are expected to be complete out to redshift z ~ 0.4 (Kim 
et al. 2001), and so should contain about 20% of the clusters in these 
deeper catalogues. 



5 X \W hr Mq'j. Below that mass threshold, the cluster pop- 
ulation expected to be detectable for Planck in the SDSS survey 
area is completely contained in the SDSS cluster sample. Since 
the cluster population is expected to die off rapidly at redshifts 
beyond unity even in low-density cosmologies, it should be pos- 
sible to identify almost all (> 98%) Planck clusters in the SDSS 
area with previously detected SDSS clusters. Thus SDSS will 
provide positions, approximate redshifts and optical luminosi- 
ties for almost all clusters Planck will see over more than a third 
of the usable high-latitude sky. Planck is expected to detect of 
order 1.5 galaxy clusters per square degree at the 2-o signifi- 
cance level. Therefore, it should be possible to identify ^ 15000 
clusters in the SDSS data which are also detectable for Planck, 
and approximately 70% of those or ^ 1 1000 clusters will be ef- 
ficient weak lenses (Bartelmann 2001). The combined data will 
allow detailed, multi-wavelength studies of a rich, uniquely and 
homogeneously selected galaxy cluster sample. Correlation of 
the two data sets (e.g. requiring detection at > 2-o in both or 
stacking Planck data at the positions of SDSS clusters) will al- 
low one to go much deeper (cf. Fig. |l^. Note that the Planck 
clusters which are not visible in SDSS should almost all be at 
z > 1.2, so that SDSS offers a way to identify the small high 
redshift tail of the Planck cluster distribution. 

4. Summary and Conclusions 

Dalcanton (1996) suggested searching for galaxy clusters in op- 
tical surveys by searching for excess surface brightness in heav- 
ily smoothed images. A recent application of this technique, the 
Las Campanas Distant Cluster Survey (Gonzalez et al. 2001) re- 
sulted in a catalog of 1073 groups and clusters with redshifts out 
to -0.8. 



Drift-scan surveys like the LCDCS are ideal for this type of 
project as they naturally yield highly uniform images. We inves- 
tigated in this paper redshift limits for cluster detection in the 
Sloan Digital Sky Survey data. We assume that 

• clusters have surface brightness profiles following the den- 
sity profile suggested by Navarro, Frenk & White (1996, 
1997); 

• the Universe is well described by a ACDM model (i2o = 
0.3, Q.A = 0.7, h — 0.7) normalised to the present number 
density of galaxy clusters; 

• the spectral energy distribution of clusters is a superposi- 
tion of a dominant early-type and a late-type galaxy spec- 
trum; 

• the early-type fraction of galaxy clusters decreases with in- 
creasing redshift, reaching half of its present value by red- 
shift unity; 

• early-type galaxies brighten with redshift; and 

• the mass-to-light ratio of galaxy clusters at present is 
25OM0/L0 in the Johnson B band. 

Then, using a measurement of the sky brightness at the 
SDSS telescope site, applying the SDSS photometric system and 
smoothing with a Gaussian kernel, we showed that: 
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• The most efficient single SDSS band for galaxy-cluster de- 
tection is the /' band, in which the 5-a detection limit for 
clusters with mass M ~ 5 x IO'^/i^'M© is approximately 
unity. 

• Summing data from different bands takes one considerably 
further. While the m' and bands are not very efficient for 
this purpose, summed r' + /' + z' data allow clusters with 
M ^ 5 X 10^'^ h^^ Mq to be detected at 5-a significance out 
to redshift ~ 1.2. 



• The limits derived here are in good agreement with extrap- 
olations from the Las Campanas Distant Cluster Survey. 

• Brightest cluster galaxies should be detectable in the SDSS 
data out to redshifts near unity. For a substantial fraction of 
the SDSS cluster sample, it will therefore be possible to de- 
rive redshift estimates from the photometry of the brightest 
cluster members associated with them. 



Finally, it is worth noting that, except for the most massive 
clusters, the SDSS cluster detection redshift limit falls above the 
upper redshift limit for cluster detections in the sub-millimetre 
regime expected with the upcoming Planck satellite. This im- 
plies that an SDSS cluster sample should contain almost all the 
clusters Planck is expected to see in the SDSS area, and that it 
will thus become possible to study a sample of more than 10"* 
galaxy clusters in both the optical and sub-millimetre regimes, 
most of which will be efficient weak gravitational lenses. 
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